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ПРЕДИСЛОВИЕ 


Цель настоящей книги — предложить инженерам-расчетчи- 
кам и студентам формулы численных квадратур (для прибли- 
жённого вычисления определённых интегралов и площадей), 
удобные для пользования и в то же время обладающие возможно 
большей точностью получаемого результата. 

Формулы, предложенные Ньютоном и вычисленные его со- 
временником Котесом, являются наиболее простыми: они имеют 
равноотстоящие ординаты и достаточно простые множители 
при них, но дают ошибку, обычно в несколько десятков раз 
ббльшую, нежели наивыгоднейшие формулы Гаусса при таком же 
числе ординат. В свою очередь, формулы Гаусса неудобны 
для пользования из-за несоизмеримых абсцисс и таких же мно- 
жителей при ординатах. 

Промежуточное место занимают формулы Чебышева, кото- 
рые имеют простые одинаковые коэффициенты при ордина- 
тах; но несоизмеримые значения абсцисс затрудняют пользование 
ими, особенно при вычислениях, требующих значительной 
точности. 

Возвращаясь к формулам Гаусса, дающим наибольшую 
возможную точность результата, приходится отметить, что вы- 
числения по этим формулам весьма громоздки, и поэтому ими 
в практике пользуются сравнительно редко. Академик А. Н. Кры- 
лов об этих формулах говорит: «При одинаковом числе орди- 
нат формула Гаусса даёт вообще ббльшую точность, нежели 
другие формулы, но зато вычисление по ней сложнее; поэтому 
её выгодно применять в том случае, когда вычисление орди- 
нат затруднительно» *). 


*) «Лекции о приближённых вычислениях», стр. 107, изд. 1933 г. 
19 » 3 
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Мы поставили задачу — дать такие формулы квадратур, 
чтрбы, не увеличивая трудоёмкости вычислений, получать pe- 
зультат возможно более точный. Чтобы достигнуть этого, мы 
берём абсциссы, близкие к гауссовым, но удобные при вычи- 
слениях, в сотых или девяностых долях интервала интегри- 
рования. Такие абсциссы дают нам и коэффициенты при 
ординатах, выражаемые в виде простых дробей, которые 
в некоторых частных случаях можно значительно упростить, 
заменяя их. приближённо близкими к ним дробями с малым 
числом. знаков и учитывая, в случае Медо врсти, весьма малые 
поправки к этим дробям. 


Ар: гст.огд.р)! 


$ 1. ОБЩИЕ СООБРАЖЕНИЯ 


Когда мы составляем интерполяционную формулу с целью 
получить из неё формулу квадратур, ординаты следует рас- 
полагать симметрично относительно середины отрезка, для ко- 
торого составляется интерполяционная формула; иначе площадь 
кривой и её зеркального отображения будет несколько различ- 
ной, за исключением случаев, когда интерполируемая кривая 
есть полином не выше (п-- 1) 
степени при л данных для фор- 
мулы квадратур. 

Замечание. Возьмём какие- 
либо две ординаты на кривой 
(черт. 1, @), расположенные по 
обе стороны от середины интер- 
вала интегрирования. Высоты ор- 
динат пусть будут аи В, абс- 
циссы их—а и 6. Если opan- Черт. 1. 
наты расположены симметрично 
относительно середины отрезка, то множители при них будут 
одинаковы; обозначим их через А. В выражение площади кри- 
вой очи войдут в виде суммы произведений (Аа-|- А8). В зер- 
кальном отображении кривой (черт. 1, 6) абсциссе а соответ- 
ствует ордината В, абсциссе 6--ордината а. В выражение 
площади они войдут в таком же виде: (А8 -|- Аа). 

Если ординат расположены несимметрично, то коэффи- 
циенты при ординатах будут различны. Пусть при ординате а, 
соответствующей абсциссе а, — коэффициент А), при opan- 
нате В,, соответствующей абсциссе 6,, — коэффициент В. В зер- 
кальном отображении абсциссе а соответствует ордината В 
с коэффициентом А, абсциссе В, — ордината а, с козффи- 


*) А, отлично от А, так как смещение одной ординаты изменяет 
коэффициенты при всех ординатах, 


2 в. П. Ветчинкин и Ф, М. Коган 5 
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циентом В. В выражение площади основной кривой эти новые 
ординаты войдут в виде (Аја ВВ), а в выражение площади 
её зеркального отображения в виде (4,8 —- Ва). Очевидно, что 
эти выражения могут оказаться равными только в очень редких 
случаях; вообще же площади кривой и её зеркального ото- 
бражения будут различны, 

Чтобы решить поставленную задачу, мы прежде всего на- 
пишем интерполяционную формулу Лагранжа в самом общем 
виде *): 

— y, 2 — Хз) (# — лу) (ж — 2)... (# — Ха) 
У Ni — xa) (X1 — X3) (X1 сл)... a — ха) + 
(x — ху) (x — ху) (x — х,)...(х — Хр) 
+ Уа = ха) (£2 — Ху) (Xa — №4)... (X2 — Хи) ет 
К-х)«--хо)«--ха)...(«--ха) 
Ера, — р) (а — Хо) (Хз — Xa) · „(Ав — Ха) 87. 
паронд (x — х\) (x — x) (£ — хз)...(* — Ха—) | 
n (£n — Х\) (Ха — Ха) (Sn — Ху)... (р р) 


Здесь в #-том члене в числителе отсутствует множитель 
(х —х,), а в знаменателе — множитель (ху — х,) =0. 

При симметричном расположении ординат и нечётном их 
числе (черт. 2) формула примет более простой вид: 


т—1 т 
G Fam Faw П (2 + а) в Писмо 
=1 == 


У = Уут + 


аан) П t am + a: аа П (— am — a) 
i= i=l 


т-2 т 

(+ ап) (а F ан). П ада. Пе сво 
i= i= 

та т 


(ат + ап) аи). [| ан, + ad-an П ата) 
1=1 = 


+ У—(т—‹) Њ 


т т-1 
Пато. в По (к о) 
+...+ж = Аы 


1 
Пал + 2) ап] • п (ат — а)-(х — ат) 
1=1 i=l 


*) А. Н. Крылов «Лекции о приближённых вычислениях», 
стр. 81, изд. 1933 г. 
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Здесь числители всюду одинаковме: 
(22 — аф) (ха — абу)... (х — а) [х] = М, 


а знаменатели у коэффициентов при симметричных ординатах 


Ф сав Фил "а 070 40 да ба h 
Черт. 2. 


отличаются лишь одним множителем, и мы можем написать: 


ју 1 № У—т т | | 
кет== НУ Скај податок Коха таб У 
7 2 | а? (а, — ай: 3) (а, = ад, да, — а [5 Ham + хлат | 
№ 
+з 


2\2 2 5 5 2 
(ар ар д-ад-1-#т--2)--44щ-1 744) 


У-(т- Ym- го а же 
хо ва | +. ааа тање Х 


Нат 0—4 m1 (а1—а5) (01—а?, 1). а 
Ул и анс аа 2% чара 
% Е се нина тт Юта? ай аа; “ ) №0 


Если п чётное (черт. 3), то в каждом члене с индексом ; 
отсутствуют соответствующие средней ординате множители [х | 


в числителе и (ај в знаменателе. Числители будут опять 
всюду одинаковые: 


(x? — ап) (x? — ат). (43 — М) = №, 4 
в | 1 
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и будем иметь: 


(си ут Бе ка заден 


Де атл Im) amla ул ЕЕЕ Ж 
У—(т—41) Ут—\ № 
x слез кк F F 1 


(а? 1—а2) (аз--а“, „4“ Ха?—а? ја, 
У У 
х [= а, Ж ар (Ш 


Приближённую величину площади интерполируемой кривой, 
даваемую интегралом интерполяционной формулы, обозначим 
через 5. Её аналитическое выражение будет иметь такой вид: 


5= Ауу, - Ау, АУ +... А„у„ (Ш) 


где все коэффициенты А, зависят только от расположения орди- 
нат у, но не от их величин. 
Чтобы иметь возможность оценить ошибку, происходящую 
от замены искомого интеграла 
+h 
F= {гах 
h 


данной кривой ero приближённым выражением 8, мы предпо- 
лагаем, что интегрируемая функция разложена в ряд Тейлора: 


х? х3 х х5 xê х1 
x10 хи хі? 
к: НЕ я бере ЕР Чес 
тогда её интеграл выражается рядом 
+һ 


| гах= роза стр? рай Яка 
идр), 


гле 
А=2 (при х<< 0; В=2; Сег D=} ит. д. 


Чтобы найти остаток R, определяемый как дополнение при- 
ближённого значения площади до точного 


Р= 5+ К, 
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разложим каждую из ординат у, (формулы (Ш)) в ряд Тейлора 
и будем сравнивать коэффициенты при одинаковых степенях И 
в выражениях для $ и для Р. 

Вообще говоря, число совпадающих коэффициентов опре- 
деляется числом ординат л и не зависит от их расположения; 
оно равно 5 или ЕЗ соответственно для л чётного и He- 
четного. Но возможно специальное расположение ординат, 
дающее совпадение ещё одного или нескольких коэффициентов 
при одинаковых степенях Å; их число определяется числом 
свободных параметров, входящих в выражения коэффициентов 
А, и зависящих от расположения ординат; число их равно 
п п s 


—1 
g или | —5— | соответственно для л четного и нечетного. 


Наивыгоднейшее расположение ординат было найдено -Гаус- 

сом, который потребовал, чтобы формула квадратур (Ш) была 
. вполне верна для всякой целой функции степени не выше 
2п— 1. При этом для кривых, разложимых в ряд Тейлора, 
автоматически удовлетворяется и другое условие, которое одно 
только нас и интересует: разность между площадями истинной 
кривой и её интерполяционного выражения (взятыми в преде- 
лах интервала интегрирования), т. е. Р — S, получается наименьшей 
из всех возможных, 

Мы подходим к разрешению интересующей нас задачи, 
рассматривая разности интегралов. Это позволяет нам не только 
отыскать условия минимума, но и найти, насколько мы от него 
отходим, отступая в ту или другую сторону от оптимального 
расположения ординат. ' 

Интегрировать интерполяционную функцию и находить оста- 
ток К-- Е--8 мы будем не в общем виде, а для частных 
значений л, так как наибольшее интересующее нас число орди- 
нат равно семи. 

Начнём с п==2. 


$ 2. ФОРМУЛЫ КВАДРАТУР ПРИ ДВУХ ОРДИНАТАХ 


Пусть абсциссы наших ординат будут: (— ай); (ай). 
Формула Лагранжа для двух ординат даёт 


у= уа (хай) у (е ай). | 
Интегрируем: 


+ в 
Цудеет зету, у) = Е), (2) 
=h | 
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Чтобы найти остаток, выражаем ординаты у, = Z; и у, == 2, при 
помощи ряда Тейлора: 


ла АЗЕ Вай СО КОИ НЕЧТО ВОТ F 
У асве вен овер.) 


подставляем в выражение для S; имеем: 
ah? athi 4886 
5==2 ater HEG 1057 Рич 
=> {arpe t 1 8а% ЕК Биро Итар...) (27) 


Вычитаем (2) из интеграла функции J, также выраженной в виде 
ряда м (1), и вх: остаток 


в--2 {с gy 01—842) E fi (1--5а9 ОА (1—7Та)-=...). 


Здесь имеется один свободный параметр а. Можно его выбрать 

таким, чтобы коэффициент при первом члене остатка обратился 
T 

в нуль; это даёт 1 — За? =0, откуда a= у — это наи- 


выгоднейшее значение а, которое соответствует формуле Гаусса, 
и для случая двух ординат — также формуле Чебышева; её 
остаток начинается с члена, и h5: 


Вместо несоизмеримого значения а== ит --0,57785... 
мы можем взять близкие соизмеримые значения абсцисс; а == 
= = 0,57777... ; 0,58; 0,6 и др. При каждом значении a 


найдём коэффициенты при ординатах и остатки. Для сравне! ия 
точности получаемых каждый раз выражений площади мы выбран- 
ные нами абсциссы, получаемые при них коэффициенты при 
ординатах и остатки будем располагать в виде таблицы. 

В случае двух ординат при. любом выборе абсцисс коэф- 


1 
фициент при ординатах ‘равны -,, остатки же пойузансо 


2 , 
различные (см. табл. 1). 
Имея возможность легко находить остатки любой формулы 
квадратур, мы“ пользуемся этим и приводим в наших таблицах 
не только формулы Гаусса и соседние с ними (т, е. при ордина- 
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тах, расположенных по соседству с наивыгоднейшими гауссовыми 
ординатами), но и все остальные формулы квадратур: Ньютона- 
Котеса, Чебышева, Понселе, Штаермана и др. 

При этом формулами типа Маклорена будем называть 
формулы, в которых равноотстоящие ординаты взяты в сере- 
дине интервалов вместо ординат на концах интервалов (черт. 4). 


-À + + 0 •ћ 
Черт. 4. Черт. 5. 


Круговыми называем такие формулы, где за абсциссы 
взяты проекции п точек деления полуокружности на равные 
части (черт. 5). 

Формулами типа Понселе назовём соответствующие 
ординатам, построенным в серединах равных промежутков и, 
кроме того, в двух концевых точках (черт. 6). 

Табл. 1 (см. в конце книги) показывает нам, что наиболее 
употребительная формула трапеций является наихудшей из всех 


Черт. 6. 


формул квадратур на 2 ординатм, не исключая даже случая 
слияния их в одну, расположенную в середине интервала инте- 
грирования, Это видно из чертежа 7. 

5 3. ФОРМУЛЫ КВАДРАТУР ПРИ ТРЁХ ОРДИНАТАХ 


Ординаты возьмём в точках (— ай); (0); ( ЄР). Здесь так- 
`же имеется один свободный параметр а. 
По формуле Лагранжа 


y = зда Wax се + ah) Руж (х — ah) — у, (ха — аа). 
и 


аи 


| дх „29% 2 ЗА 25 tay Фау, | = 


= (роти 2 (#—)») - (8) 


Разлагаем ординаты в ряды Тейлора (причём здесь и далее 


м 


аа. ами 


13 == А Вал 4- Ст 2057 Е РФ и | 
За 1 848 — 
рое НТК... 
у: =A, 
4 
иа (АС вет ока" +), 
2y, 2А. 
Подставляя (Y-Y) и Уз в сре для зае имеем: 
2 
= {2 (а ст ДЕСИ +09 ком +...) + 


ао па) (асвет за y gt Tet 
Ка +..). (8) 


Вычитая (3’) из интеграла функции г, тоже разложенной в ряд 
Тейлора (1), получаем остаток 


в- Е (1—) +07(1—%) + Ка (1—%)-+...} 


Чтобы обратить в нуль коэффициент при первом члене остатка 
1— за ‚ нужно положить а == 5 ; тогда получаем формулу 
Гаусса; ев остаток 
т ћУ 
R = 0,32G zr t 9704К gF .. 
Для формулы Чебышева должно выполняться равенство козф- 


фициентов при ординатах ==? (8-3); откуда а== FT А 
Остальные формулы не требуют пояснений; относящиеся к ним 
коэффициенты приведены в табл. 2. 

Наиболее близкие к гауссовым значения абсцисс a=?!) o 
дают наименьший остаток; но абсциссы 1), выгоднее тем, что 
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коэффициенты при ординатах в этом случае более простые, 


остаток же увеличивается незначительно, 


Из табл. 2 мы снова видим, что из всех формул квадратур 
на три ординаты наиболее употребительная — Симпсона — по 


величине остатка является наименее точной. 


$ 4. ФОРМУЛЫ КВАДРАТУР ПРИ ЧЕТЫРЁХ ОРДИНАТАХ 


За абсциссы примем (— Ёл); (— ай); (ай); (bh). 


Интерполяционная формула для четырёх ординат: 


р пива = щие Ами bh) — y; (x — bh)] — 
— b (x? — 5282) [ у, (хай) — у, (х — ап)]). 
Интегрируем: 
+h 
| y dx =S = 


= зит (Ф аи 1) а уд] И 


Поступал аналогично предьдущему, находим остаток; 


ЗИД Өр дилање pa 


on [| — Тци-ри-ав) — зар (+) 4 


+x и gatat ане) 


—3a?b? (at 4 bt 4 а) |. > 


Здесь имеем два свободных параметра: а и b; следовательно, 
можно избавиться от двух первых членов остатка. Чтобы обра- 


тить в нуль первый член, нужно выполнить условие 
1—2 (a pb? — 30252) = 0. 


Отсюда находим соотношение между а и 6: 


3 — 582 
а= ОНО 
ет! 
Џи. 2582 . 

5— 154 


Чтобы (тен и второй член, потребуем, чтобы 


тари 4 аза) — Зайра (a? 40) == 0. 
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(5) 


(57 
(5) 


(6) 
13 


Из (5“) и (6) получаем уравнение 
5,654 — 4852 --0, 48 ==0, 


откуда находим: а==0,339981; 5 = 0,861136 — это абсциссы 
формулы Гаусса. 

Вместо несоизмеримых значений абсцисс Гаусса мы берём 
соизмеримые значения, представляющие или двузначную десятич- 
ную дробь или простую дробь с малым числом цифр в числи- 
теле и знаменателе. Дроби со знаменателем 90 особенно удобны 
для вычисления площадей функций, имеющих аргументом гра- 
дусы окружности. Задавая значение одной из абсцисс, мы нахо- 
дим по формуле-(5' или (5") наивыгоднейшее значение второй 
абсциссы несоизмеримым, после чего подыскиваем ближайшее 
к нему соизмеримое, на котором и останавливаемся. 

Наиболее близкие к гауссовым значения абсцисс 


а—0,34; 6--0,86; 


или 
а= {уу==0, 338888. b= S=, 861111. 
Им соответствуют множители при ординатах: 
а = 0,3255342 и =, 1744658 в первом случае 
и 
13 225 7079 
40608 = 0,325674 и 40608 =0,174326 — во втором. 


В том и другом случаях множители при ординатах с очень малой 
неточностью можно заменить простыми дробями: 


28 > 5. 
gg == 032558 14 и 86 == 0,1744186, 
а также ещё более простыми; 
13 7 
40 = 0,325 и 46 = 0,175. 


Чтобы не вносить этим ошибки, можно к упрощённому множи- 
телю добавить малые поправки, данные в табл. 3. Формулы 
тогда примут вид: 


5=>2- (5 (у; Љу Ри +у)— 
— 0,0000933 [(у; + ya) — (у, + У] ) , 

5 —=2/Л (5 у) До +) + 
+ 0,0006747 [(у, -- и) — (у; +») + 


Малыми поправками можно пренебрегать, если это допускает 
желаемая точность результата. у 

Кроме упомянутых здесь наилучших приближений, мы даём 
в табл. 3 также более далёкие от наилучших, но проще выра- 
жаемые абсциссы и множители при ординатах. 

Для формулы Чебышева должно быть 


откуда 


Вторым соотношением берём (5), которое с предыдущим равен- 
ством даёт: 


2 2 де 
3—5 [e+ (+—е) — за (3 а) | =, 
Получаем уравнение 
4544 — 30а? -- 1 ==0, 
из которого находим 
а==0,187592 и ђ==0,794654, 


т. е. приходим к абсциссам формулы Чебышева. 

Вместо полученных несоизмеримых значений берём соизме- 
римые. Положив b= 0,80, находим, что a= 0,16330 (если одна 
из абсцисс удаляется вправо, то вторая уйдёт влево). Чтобы 
ещё более упростить формулы, мы можем взять для а значе- 
ния 0,16 или 14. Остатки, приведённые в табл. З, показывают, 
насколько портится результат от таких упрощений. 

Профессор И. Х. Штаерман сделал попытку приблизить 
формулы Чебышева к формулам Гаусса. Для этого он равные 
множители при ординатах Чебышева заменил однозначными мно- 
жителями, близкими к таковым же (но несоизмеримым) в фор- 
мулах Гаусса. Абсциссы при этом, как и в формулах Чебышева, 
получаются несоизмеримыми. Например, в случае 4 ординат 


1 
Штаерман заменил множители 40 и 40—40 Хт ль со 


Ти 2=78 . Тогда из общей формулы (4) находим 
1 1 
:(:—е)=(н– 1), 
откуда 
г 2—1 — 24а, 
или 
2а? 4 2—1. 
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Полагая злесь 4-- 0 и замечая, что обе средние ординаты у; 
и у; сливаются в одну, приходим к формуле Симпсона. 


Полагая а= 5, найдём 2 = 0 и, замечая, что обе орди- 
наты с абсциссами — Ó сливаются в одну, приходим к формуле 
Чебышева на 3 ординаты. 


Полагая a= b = У +. приходим к формуле Чебышева на 
две ординаты. њу 
Вообше, задавая для а любое “значение << У 5 ‚ найдём 


отвечающее ему значение 6 и величину остатка. Чтобы найти 
наивъгоднейшее соотношение между а и b, обратимся к усло- 
вию (5) и из него найдём: 


3— 5[а#-{- 1 — 2a? — да? (1 — да?) =0. 
Это даёт нам уравнение 
15а% — 1042 + 1 ==0, 


корни которого a= 0,350022 и b= 0,868890 совпалают. с абс- 
циссами формулы Штаермана. Но мы можем взять значение 
га =0,35, тогда b уйдёт вправо и будет равно 0,868907 = 
== 0,869 = 0,87. 

Рассмотрим случай, когда две ординаты построены в кон- 
цевых точках интервала интегрирования; этому соответствует 
значение параметра 6, равное единице. Имея теперь один сво- 
бодный параметр а, мы можем его выбрать таким, чтобы коэф- 
фициент при первом члене остатка стал равен нулю. Подставив 
в (5) значение b= 1, будем иметь 


1—2 (2: +1 —832:)==0 


и нахолим, что а= У 0,2==0.447218. Это даёт формулу, ne- 
давно предложенную Я. С. Безиковичем *), полученную им из со- 
вершенно других соображений. 

Кроме этого наивыгоднейшего, но несоизмеримого значения, 
мы берём близкие к нему соизмеримые; а-- 4; а==0,45; 
а==0,5. Получающиеся множители при ординатах и остатки 
помещены в табл. 8. 


+) Трулы Ленинградского индустриального ин-та, № 4, 1937 r. 
вый. П, стр. 12. 
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$ 5. ФОРМУЛЫ КВАДРАТУР ПРИ ПЯТИ ОРДИНАТАХ 
Ординаты расположим в точках: 
(— 28); (— ай); (0); (4 ал); (+ bh). 
По интерполяционной формуле Лагранжа имеем: 
Ури (ИРА (х? очна a?h?) рис 4- bh) + 
+ уе— 20) — benta (x2 — ват) у, (x+ ah) + у, (e—a) -+ 
+ 213 (b? — а?) у, (х2 — ah?) (x? — бара)! : 


Интегрируем: 
и: 2A 
>h 


X (у у) + 2 (2 — a?) [3 — 5 (a + 0°) + 1542041 уа). (7) 
Как и в предьдущем, получаем остаток: 
в=2{о| 1— двери) — 50) | az 
ка [1— i {3 (at Ен а) — барат | + 
+M in| i—i (810° -H0 еве 4-09) — 
— ба (at -+ bt -+ ал) +9 


Чтобы получить абсциссы Гаусса, коэффициенты двух пер- 
вых членов остатка приравняем нулю: 


1 — 513 (а" +0?) — 5096] =0, (8) 
1— 56804-08 + atb?) — ба (а? ва) ==0. (9) 


Из (8) находим 


МЕТР ~ 
или 

wai Веће ` 

0870603) (8") 


Соединяя (8’) и (9), получаем уравнение 6354 — 7052 -+ 15 = 0; 
его корни 6==0,906180, а==0,588469 совпадают с абс- 
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циссами Гаусса. Вместо этих несоизмеримых значений абсцисе 
мы берём близкие к ним соизмеримне; 


97 163 47, 
a= јо b= 150 или а== о; 


или a= 0,52; ђ==0,9 и т. п. 
Для формулы Чебышева должно выполняться равенство коэф- 
фициентов при ординатах: 


а? (8 — 5а?) =b? (502 — 3) = 2002 — a?) [3 — 5(b? -+ a?) 4 154222]. 


Эти равенства дают нам два уравнения для определения а и 2: 


81 
0== бо 


а (3 — 5а?) = b? (502 — 3), (10) 
62 (562 — 3) == 2 (b? — a?) [3 — 5 (5: 4 a?) + 154202). (11) 

Из (10) 
Е ВН ИШЕ Ја ат (10°) 


(10') и (11) сводятся к уравнению 
Е там — 6002 + 7 —0, 
решая которое, находим: 
b = 0,8324975 и a= 0,374541 


— абсцис̧сы формулы Чебышева. 
Пусть две ординаты построены в концевых точках, тогда b =1. 
Приравняв нулю коэффициент при первом члене остатка: 


1— 5[8 (2 1) — 59] =, 


найдём наивъгоднейшее значение параметра а, равное ИЗ 


= 0,654654, что даёт вторую формулу Я. С. Безиковича *). Наибо- 
лее близкие к полученному соизмеримне значения а= или 


а —=0,65 дают малую ошибку; но мы берём и близкие чётные 
значения а==0,64; а==0,66. Соответствующие множители при 
ординатах и остатки даны в табл. 4. И. Х. Штаерман для 5 
ординат даёт формулу, где множителями при ординатах служат 
дроби: 

3. 4 469-300 3 

пт“ очи 
и получает, что значения абсцисс должны быть а == 0,546377; 
р = 0,905335 (в значение крайних абсцисс у автора вкралась 


Ы Труды Ленинградского индустриального ин-та, № 4, 1937 г,, 
вып, П, стр. 13. 
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ошибка; точное вычисление при так заданных множителях даёт 
b = 0,9053237). 
Мы предлагаем аналогичную формулу с множителями при 
ординатах: 
{- Ий. жылу > 
25" 25' 25’ 25’ 25" 


Ординаты при этом берутся в точках: &= 0,5337205; 
b= 0,9050822. Такая формула даёт лучший результат (мы далее 
приводим пример), а главное, удобнее оперировать с дробями 
в двадцатьпятых долях. 


$ 6. ФОРМУЛЫ КВАДРАТУР ПРИ ШЕСТИ ОРДИНАТАХ 
Пусть абсциссы будут: 
(— ch); (— bh); (--ай); (Дай); (-- 51); (ей). 
Интерполяционная формула для шести ординат: 


она е Да АА шл. му 
У == habe (с — ай) (с — 52) (62 — a3) 


X {abh (b? — a?) (x? — ат) (x3 — b?h?) уж Ней) — у (е — 
— сћ)) — ach (c? — а?) (x? — а?) (x? — eh?) [y(x + bh) — 
--у(х--0й) = бей (c? — b?) (x? — 22) (х2 — с?Л®) |у(х-ай)-- 
- —(®* — ah))}. 

Интегрирусм: 


+, 
ра раа Х 
х (ао Ну © + ав) - 
ео нь |526] + 
+05 "+ на || = 


ња ЕВ (Сл + у) (02 — а*[3—5 (а 4-62) 4- 


+ 159262] — (у, -|- уь) (c? — а) [3 — 5 (a2 -H e) 4- 15а2с2]4- 
Јоце в) [3—5 (62 4- се) 1592 ]}. (12) 
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Поступая так же, как и ранее, получаем остаток: 
= 200 9011—8030 00) Бато аа 4 ма) + 
4 1292) | 4 
+ В {3 (at 4 bt 4 ct аз ае A ре») — 
— 5 (a? $ b?) (a2 4 с^) (b? 4 c?) лага (а 4 з +e] + 
+ [1—15 BeH e ato ра) аш] — 
о Ка“ -+ аќс& + біс“) + (а262 + а?с? + Р2с2)(а* + 4- + 
-+ 2262с2(а? -|- 62 + с? 1 + 1562224 54 — сї -p 4202 -|- ас? -- 
+ ма) + Р 1—13 (3 [abt ал не 
фка 4 раја езу] — 
— 5 (42 b?) (a? + c?) (024 с2) (at -H01 -+ ct) 4 
-H 15a2b2c? [(а 4 м 4c?) (at bt ct) 4 ара) | +... 


Теперь у нас три свободных параметра: 4, 5 и с, что поз- 
воляет избавиться от трёх первых членов остатка. Приравняв 
нулю коэффициенты при них, мы получим систему трёх урав- 
нений, решение которой даст нам наивыгоднейшие значения 
а, би с— абсциссы Гаусса. Но мы проделывать эти громоздкие 
для шести ординат (и аналогичные для семи ординат) выкладки 
не будем: нет смысла искать уже известные абсписсы, тем более, 
что у Гаусса это делается проще — им найдена общая формула, 
позволяющая найти абсциссы при любом числе ординат. 

Это же можно сказать и о формуле Чебышева, абсциссы 
которой мы нашли бы, приравняв коэффициенты при ординатах. 

Здесь мы ограничимся отысканием наивыгоднейших значений 
двух параметров, считая третий заданным. 

Пусть c= 1, т. е, две ординаты построены в концах OT- 
резка. Тогда мы можем два первых члена остатка приравнять 
нулю: 


1— ЕС -+ b? + с2)— 59262 -4 a?e? |-62с2) -+ 15а262с?] = 0, (13) 
j 1—5, [3 (at bt- ct- a?b? ае ма) — 


—5(a? 4 (а? суб -+ с) 4 15a2b2c? (а? +b? + с?у] == 0. (14) 
29 


Эти уравнения, после подстановки C=] и исключения а, CRO- 
дятся к бикубическому уравнению 


10568 — 91654 | 1562 — 1 ==0, (15) 


которое заключает лишний корень, соответствующий случаю по- 
парного совпадения ординат у, су, и у, с Уу, что приводит 
нас к случаю четырёх ординат, из которых две расположены на 
концах отрезка. Мы уже видели, что средние абсциссы тогда 
равны + У 0,2. Разделив уравнение (15) на (b? — 0,2), nony- 
чим биквадратное уравнение 
10554 — 7062 + 5 =0, 
из которого находим 
0--0,285232 и 6==0,7650552. 


Но мы округлим полученные числа до сотых долей и положим . 
a= 0,28; b = 0,76. При этих более простых абсциссах первые два 
члена остатка не исчезают, но коэффициенты при них получаются 
малыми; вычисления же при таких абсциссах много проще. 

Если мы зададим два параметра, то можем найти наивыгод- 
нейшее значение третьего. Эти два параметра положим близкими 
к гауссовым значениям абсцисс — пусть 


а=0,24; 6-- 0,66. 


Имея теперь один свободный параметр, мы выбираем его таким, 
чтобы коэффициент первого члена остатка обратился в нуль. 


Из (13) получаем: 
са — 15—21 (a 463) 4-35 ама 
21 — 35 (а 52) + 105а262“ 
Подставив сюда значения 
а=0,24; 6--0,66, 
найдём, что 
с? —=0,866373005 и c= 0,9307916. 
Мы берём значение с = 0,93 и, кроме того, чётные: с = 0,932 и 
е == 0,94. Соответствующие коэффициенты при ординатах и остатки 
помещены в табл. 5. 
Формула Штаермана для 6 ординат имеет множители при 
ординатах: 
а Вада. 
12; В; боб 
а абсциссы расположены в точках: 
a= 0,159165; b= 0,681582; с == 0,512038. 
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$ 7. ФОРМУЛЫ КВАДРАТУР ПРИ СЕМИ ОРДИНАТАХ 
Имеем семь ординат с абсциссами 
(— ch); (— БА); (— ал); (0); (Чай); ( -+ bh); (+ ch) 
Интерполяционная формула 


1 А 
V= уа) Х 
х (x? — а) (x? — 6282) [у, (х + ch) + y, (х — ет) — 


x (а — а? 
— © x (x? — ah?) (x° — сив) [ye (4-20) + ya (ес) + 
x (x? — бен) (8 — от [у, (x ед) -- у(х — а) — 


4 aze 
„(е — а?) (c? — b?) (b? — а?) у, (x? — ah?) (x? — боја) о) 
Интегрируем: 
+ : 
2 — а? 2 2 
уе (т S Ttt + D)o Hy 
Е 9 
(Ее) 
(сї — а?) (c? — b?) (b? — аз) (1 сыны 4 
айг одо 
4. таа КАЖИ ара) А 2v} = 
= аара r је с? “15—21 (а? -+ 52) + 
бану, у) — 56015 — 2162-4-05) ване + у 
ŽE 15-01000 4-е) зам) (у, Ву) — 
Е а (е ие) 
+ 35 (4262 А- a?e? |-62с?) — 105а20:с2]. ду, ) . (16) 
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Поступая аналогично предыдущему, получаем остаток: 
в-4ка |! = КРИТ са) — 21 (ара ага ра). 
-H35a2b2c2} | ци! iii ЧЕ са 1544-9444? јазел. ф2е)-- 
— 21 (a? 4 b?) (a? 4 3) (b2 -4 c?) 4 35a?b2c3 (a? ија) ] 18 
Рав (19064-08 A atot et) а] — 
соја ас 4 btet) (и је 4 bac?) (at 4 bt ед. 
4 абс? (a? 4 b? 4- с2)] -4 35a?b2c? [at A bt -p ct 4 азђа аа | 
+ ое) | леа! — 19 (15Щачи ае меч 
е-рацаца-рица-раие)-- 

о-в) (еее 

4- збазраса ја 4b? 4 се) (at 4 bt ct) ова | АМА 


Приравняв нулю коэффициенты двух первых членов остатка, 
найдём уравнение для отыскания наивыгоднейших значений двух 
параметров из трёх: 


1 — 85 (a? + b? -|- сї) — 21 (4202 -|- а?с | ре) 
~ 35а242с?] =0. (17) 
1-- ыз [15 (at + bt -4 ct аго2 -|- a?e? — фе) — 
— 21 (a2 в) (a? -+ c?) (b? 4- e?) 
—- 35а202с2 (a? +b? + с2)] == 0. (18) 
Положив здесь с==1 и исключив а, получаем уравнение OTHO- 


сительно ф: 
23156 — 3096% -|- 12562 — 15 = 0. (19) 


Это бикубическое уравнение, как и в предыдущем, содержит 
лишний корень, соответствующий случаю попарного совпадения 
ординат у, с у, и уу с Yg что приводит нас к уже рассмотрен- 
- ному случаю пяти ординат с ординатами на концах, когда абс- 
циссы ординат у, и у, есть х+у 3. Разделив уравнение (19) 
на (02---8/.), получим биквадратное уравнение 


231654 — 21002 -- 85 ==0, 


из которого находим b= 0,83022; q= 0,46885. 
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Мы берём значения c= 1; == 0,83; а==0,47 и 
си 1; b=" ве ро 


при которых вычисления проще и остатки достаточно малы. 
Теперь зададим два параметра и из условия (17) найдём 
наивыгоднейшее значение третьего: 
pa 35 — 45 (a? е) 4рбзаза 
45 — 63 (a2 -F са) + 10522" 
Значения двух параметров возьмём близкими к гауссовым абс- 
циссам, которые равны 0,40584516 и 0,94910792. Мы положим 
а=0,4; с==0,94 и а==0,4; с==0,95. В первом случае най- 
дем, что b= 0,719475, и возьмём b= 72; во втором находим — 
—0,740502 и, кроме —0,74, берём ġ= 0,75, чтобы 
иметь возможность делить интервал интеграции на двадцать 
равных частей. Соответствующие коэффициенты при ординатах и 
остатки см. в табл. 6. 


$ 8. ПРИМЕР ПРИМЕНЕНИЯ ФОРМУЛ КВАДРАТУР 


В качестве примера вычислим площадь > полуокружности 
Кх) =у=ү 1— x? (см. черт. 12 на стр. 58), воспользовавшись 
формулой типа Гаусса на 5 ординат с абсциссами @= 0,52; 
= 0,9 (формула 7, табл. 4)“). 

Формула запишется: 


S= 2h [0,1256837 (у, + у;) + 0,2898774 (у, + у) + 
4-0,2688777у,), 


причём в нашем случае k= 1. 
Кривая наша симметрична, что уменьшает вычисления, так как 


У u кууу, 
yı = ys = /(0,9)= УТ—0,81= И 0,19 = 0,4358899, 
уо у, = f (0,52) =У1— 0,2704 = У 0,7296 = 0,8541662, 
Уз =f (0) 1. 
Подставив полученные величины ординат и #=1 в формулу, 
имеем: 


5 =2(0,1256837 · 2. 0,4358899 -+ 0,2398774- 2- 0,8541662 + 
-} 0,2688777) = 2- 0,788237 = 1,576474. 


*) Если бы надо было вычислить площадь какой-либо кривой, Bbl- 
ражаемой тригонометрической функцией, то улобнее было бы восполь- 
зоваться формулой, где абсииссы ланы соизмеримыми с градусным 
делением окружности, например формулой 8, табл. 4, 
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4+1 
Точное вычисление интеграла | ИТ-- х? 4х даёт 1,57080, 
= 


так что ошибка равна 0,00567. 

Посмотрим, как изменится результат, если мы возьмём при- 
ближённые значения коэффициентов при ординатах (формула 7, 
табл. 4); 


$=2(0,125.2.0,43589 -+ 0,24.2-0,854166 + 0,27) = 
=2.0,788972 = 1,577944. 


Ошибка теперь равна 0,00714, т.е. увеличилась в 1,3 раза. 

Заметим, что так как расхождение точного результата и при- 
ближённого начинается уже в третьем знаке, то можно было 
брать ординаты с четырёхзначной точностью. 


$ 9. СРАВНЕНИЕ ТОЧНОСТИ РАЗЛИЧНЫХ ФОРМУЛ 
КВАДРАТУР 


Оценка точности формул квадратур по остаткам является и 
затруднительной, так как требует знания производных высоких 
порядков от интегрируемой функции, и во многих случаях не- 
возможной, когда эти производные принимают нулевые или бес- 
конечные значения. Поэтому мы в дополнение к таблицам остат- 
ков проведём сравнение формул квадратур на ряде специально 
подобранных примеров. 

Простая синусоида у = sinx не представляет никаких OCO- 
бенностей. 

у == ѕіп2х и у = sin? х *) имеют весьма малые ординатн у KOH- 
цов кривой и довольно крутые уклоны в промежуточных точках. 

Полуокружность у= У 1— х? имеет на концах вертикальные 
касательные и соответственно бесконечные значения производных. 


Кривая у= ут имеет на концах ординаты, стремящиеся 


к бесконечности; интеграл этой кривой несобственный, но имеет 
конечное значение. Все формулы с ординатами на конпах к Bbl- 
числению этого интеграла непригодны, тогда как формулы типа 
Гаусса и Чебышева дают результаты хоть и заниженные, но 
довольно близкие к истинному. 

y=|xV x (1— х) представляет несимметричную кривую, 
встречающуюся при вычислении моментов инерции самолётных 
стоек. 


*) Аналогичные кривые встречаются при вычислении моментов 
инерции самолётных стоек и винтовых дужек. 


25 


http://rcin.org.pl 


Подинтегральнье функции зллиптических интегралов первого 
и второго рода: 


прера " умне 


при модуле #, близком к единице, являются образцом несимме- 
тричных кривых, имеющих весьма большую кривизну при х, 


7 
близком к 5. 


Мы приводим также пример Штаермана, так как в нём произ- 
ведены полные вычисления по формулам Гаусса и Чебышева, 
что мы в своих примерах делаем лишь в тех немногих случаях, 
где вычисления не являются громоздкими. 

Формулы с неравноотстоящими ординатами весьма хороши 
при вычислениях определённых интегралов, но мало пригодны 
для вычисления интегралов неопределённых и, особенно, для 
численного интегрирования дифференциальных уравнений. В этих 
случаях формулы Ньютона-Котеса являются незаменимыми. 

Для каждой из взятых нами кривых мы составляем таблицы 
площадей, вычисленных по различным формулам квадратур. 
Кроме того, для каждой кривой мы формулы располагаем в таб- 
лицы в порядке уменьшения даваемой ими ошибки. Мы видим, 
что такое расположение формул для каждой кривой своё; но 
почти всегда формулы типа Гаусса дают меньшую ошибку, чем 
другие формулы квадратур при том же и даже при большем 
числе ординат. 

В отдельных случаях менее точная формула даёт лучший 
результат, чем более точная. Так, например, для эллиптиче- 
ского интеграла | рода формула Симпсона, повторённая два раза 
(5 ординат), дала лучший результат, чем более точная формула 
Гаусса на 6 ординат. Подобные результаты объясняются 
случайной компенсацией ошибок. 

Иногда приближённая формула может дать точную площадь 
кривой. Для наглядности рассмотрим площадь кривой (черт. 8), 
вычисленную по формулам квадратур на 2 ординаты. 

Формула трапеций даёт площадь кривой с избытком. Пусть 
точки Аи В приближаются по кривой от её концов к середине — 
точке 0. Формула средней ординаты (двух слившихся) даёт площадь 
кривой с недостатком. По непрерывности на кривой найдутся 
такие точки А' и В’, что, подставив их ординаты в формулу 
квадратур, мы получим точную площадь кривой. Но эти точки 
А" и В" будут совпадать с точками Гаусса только в том слу- 
чае, если изображённая на чертеже кривая будет алгебранче- 
ской не выще 8-й степени, ; 
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Заметим, что для уверенности и правильности полученной по 
формуле квадратур площади нужно знать, хотя бы примерно, ход 
кривой. Если ход кривой неизвестен, то для гарантии правильности 
полученного результата площадь нужно вычислять по несколь- 
ким различным формулам, таккак при 
пользовании только одной формулой мо- 
жет оказаться, что входящие в неё орди- 
наты не характерны для данной’ кривой. 
Так, академик Крылов строит кривую *), 
для которой формула Гаусса на 5 орди- 
нат даёт площадь, равную нулю, тогда 
как в действительности она положительна. 

Заметим, далее, что при вычисле- 
нии площадей несимметричных кривых, по- Черт. 8. 
добных кривой у = [xV ха =), мо- 
жно интервал интегрировании разбивать на два и отдельно въчис- 
лять площадь каждой части; так, для указанной кривой въчислить 
площадь от нула до 0,4 или 0,5 по формуле с малым числом 
ординат и остальную часть кривой по более точной формуле. 
Это значительно увеличило бы точность полученного результата. 

Сравнительные результаты вычислений квадратур для указан- 
ных выше примеров функций сведены втабл. 7—15. 


ТАБЛИЦЫ 


Примечание к таблицам формул 

Формулы квадратур расположены в таблицах по числу ор- 
динат. В таблицах даны положения ординат на отрезке 
[— 5: -- 1] и на отрезке (0; 1] и соответствующие коэффициенты 
при ординатах; кроме того, в правых столбцах таблиц даны 
коэффициенты трёх первых членов остатка. 

Коэффициенты при ординатах помещены под абсциссами на 
отрезке [0; 1]. Для удобства пользования коэффициенты иногда 
даются и в виде простой дроби и в виде десятичной; в этих 
случаях для краткости записи мы коэффициенты простыми дро- 
бями пишем у левых ординат, а равные коэффициенты при 
симметричных им правых ординатах записываем десятичной дро- 
бью. Когда мы даём и приближённые значения коэффициентов 
при ординатах, то располагаем их непосредственно следующей 
строкой под точными значениями этих коэффициентов, без 
повторения строки значений абсцисс. 


*) «Лекции о приближённых вычислениях», стр. 119, изд. 1933 г. 
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Четыре 


Абсциссы на 
: я Абсциссы на отрезке 
= отровно коэффициенты 
а b 
1| Ф-ла Чебы- | 0,187592 0,79 0, 102673 0,406204 
шева 1/4 1/4 


| Ф-ла типа 
Чебышева 


9200/35754 


0,330009 
0.326073 


4 | Ф-ла Гаусса | 0,339981 0,8611363 


0,0694318 
0,173997 


5 | Ф-ла типа 


Гаусса 0,34 0,86 0,07 0,33 
1633/9360= 3047/9360 = 
5'| То же ==15/86+ =28/86— 
+ 0 ,0000472= — 0,0000472= 
| › =7/40— =13/40+ 
—0 „000581395 + 0,000581395 
› 61/180 155/180 | 25/360 119/360 
7079/40608= | 13225/40608= 
= 15/86— =28/86- 
—0,0000933— +0,,0000933— 
=7/40— = 13/40- 
—1/1482 + 1/1482 


*) См. сноску на стр. 34. 
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Таблица 3 


ординаты 
Остаток: 


коэф. при Бра. арн коэф. Я 


[0; 1] и соответствующие 


при ординатах hë 
Ев бл ка & 
0, 17879 0, и 
0 0,237 0,608 
0, 19460 0,9 
1/4 0,0516 0,164 0,489 
му 0,9 
1/4 0,0516 0,164 0,489 
Пее РАДИ НИ 0,9 
1/4 1/4 —0,0259 0,0824 0,245 
0,58 0,9 
1/4 1/4 —0,0256 0,0825 0,245 
0,9 
0, Це 0,24444 --0,0107 0,205 0,524 
0,605 0,9 
0,257313866 0,242686133 0,00191 0,218 0,534 
0,669991 0,9305682 
0,326073 0,173927 0 0 0 
0,67 0,93 
0,325534189= 0,174465811 = 0,00431 0,00959 0,119 
=0,325581395— ==(0,174418605 + 0,00241 0,00506 0,0595 
—0, 000047 2= + 0, 0000472== 
=0,325+ =0,175— —0,00695 0,00178 —0,153 
+ 0,000581395 — 0 000581395 
‚ 241/360 335/360 i 
0,325674744= 0, 174325256 = —0,00293 —0,00391 0,101 
ом =15/86— —0,00197 || —0,00248 0,0498 
0, 0000933 = — 0,0000933= 
29 „325+ =0,175— —0,00509 --0,00579 0,0467 
+ 0,0006747 — 0,0006747 
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Четыре 


Абсциссы на 
; > Абсциссъ Ha отрезке 
5 отрезке 1—1; +1] коэффициенты 
2 а b 
7 | Ф-ла типа 1/3 77, 13/180 1/3 
Гаусса РУ 1800/10058 3229/10058 
8 | То же 1/3 78/90= | 6/90= 1/15 1 
| | = 13/15 25/144 (и 
9 › 31/90 78/90 12/180 59/180 
1739/10246 3384/10246 
10 > 0,4 0,9 0,05 0,3 
4/30 11/30 
11 Ф-ла Штаер- | 0,350022 0,868890 0,065555 0,324989 
мана 1/6 2/6 
11'| #) -na типа 0,35 0,87 0,065 0,325 
Штаермана 1/6 2/6 
119 То же 48/137 119/137 18/274 89/274 
1/6 2/6 


12 | Ф-ла Макло- 


Круговая 0,3090170 = 0,8090170= 0,0954915 


=-с0$ 72° =cos36° | 0,212732195 


Ф-ла Котеса 1/3 


9) см. сноску на стр. 34. 
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ординаты Продолжение 


Остаток 
0; 1] и соответствующие J 
ри ората у Ши ме” коэф. при 
5T п Кат 
2/3 167/180 
0,321037979 0,178962021 0,00302 0,0791 0,149 
2/3 84/90= 14/15 
47/144 25/144 —0,0395 —0,0728 0,00893 
121/180 168/180 
0,330275229 0,169724771 —0,00804 | —0,0292 0,00491 
0,7 0,95 
0,366667 0,133333 0,06266 —0,0261 —0,0749 


0,675011 0, во 
2/6 —0,0253 0,0500 


6545085 0,9045085 
0,212732195 
1 
3/8 1/8 —0,593 


http://rcin.org.pl 


Четыре 


= Абсциссы на Абсциссы на отрезке 
5 отрезке [—1; +1) коэффициенты 
2 а b 
1 Наивыгоднейшая| У 0,2= 1 0 | 0,2763932 
при б=1 —0,44721360 1/12 5/12 
Близкая к наи- 0,45 1 0 0,275 
выгоднейшей 157/1914= 800/1914= 
при b =1 =64/780— =326/780-- 
16») То же —0 „0000241 +0 „0000241 
> 4/9 1 0 5/18=0,2777 
11/130 54/130 
> 0,4 1 0 0, 
13/126 50/126 
19 | Ф-ла типа Пон- 0,5 1 | 0,25 
селе **), она 1/18 8/18 
же круговая 
при 6=1 | 


*) В остатках формул, отмеченных этой сноской (см. прим. к табл. 2), 
у формулы 1” —0,00167 у 
> > 11 


У 0,00160 

> > 5' 0,000177 

Е: 
формулы 6’ я 

Aae бе. | "ова 

> > и” —0,00190 

> > 16" —0,00115 


**) Получается из формулы Котеса на 5 ординат, если пропустить 
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орлинаты Окончание 
Остаток 


[0; 1] и соответствующие у я 
соня соб a Ка, при кат при 


0,7236068 1 
325/780 65/780 


1 . 
0,0820272 = 0,0167 —0,394 —0,978 
=0,0820513— 


0,0000 | 0,00810 | —0,197 | —0,510 


—0,0165 "0,459 1,069 


- | 
0,103175 а на —1 724 


при су будут следующие коэффициенты: 


среднюю ординату. 
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Абсциссъ на отрезке 
коэффициенты 


Абсциссы на отрезке 
+] 


О 


а b 


Ф-ла He- 0,832498 


бъппева 


0,374541 0,083751256 | 0,312730 
1/5 1/5 


Т | 


2 |Ф-ла Мак- 0,4 0,1 0,3 
275/1152 100/1152 


0,5=cos 60° V= 
=00,8660254 = 
=<08 30° 


3 Круговая 


0,0669873 
14/90 


4 |Ф-ла типа 


=96/700— 
--0,0004554 


= 164/700+ 
+0,0003810 


4'| *) То же 


1/18 1/4 
2187/16000=| 176/750= 


0,055 
7/51,5276892 


7 > 0,52 0,9 0,05 0,24 
х 0,125683710=| 0,239877446= 


T ју? =0,125+ =0,24— 
| +0,000683710 | —0,000122554 


8 › 47/90= 81/90=0,9 
=0,5222 | 


*) См. сноску на стр. 40, 
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Таблица 4 


Остаток 


ГО; 1] и соответствующие 
при ординатах 


0,687270 | 0,916248743 
1/5 1/5 


0,7 0,9 
100/1152 | 275/1152 | 0,145 | 0,572 | 0,864 


0,75 0,9330127 
18/90 14/90 0,0750 | 0,245 | —0,569 


3/4 17/18 
0,2346667= | 0,1366875= | 0.00947 | 0.0492 


рея =0,2342857+ |=0,13714 0,00167 | 0,0214 
00014 88 = + 0,0003810 — 0,0004554 

= °2571429-+ 

+ 0,0001488 


0,945 
0,135349290 | —0,00624| 0,0416 


0,5 0,76 0,95 
0,268877688= | 0,239877446= 0, 125683710= |—0,00308| 0,00614 
4 =0,27— =0,24— =0,125-- K 
— 0,001122312 — 0,000122554 | -+ 0,000683710 


1/2 137/180 71/180 
145787/536787=| 0,238865528 0, 125638523 0,000738 | 0,0101 | 0,0512 
=0,271591898 


4 В, П, Ветчинкин и Ф, М. Коган 37 


http://rcin.org.pl 


Пять 


= Абсциссы на отрезке Абсциссы на отрезке 
+ [—1; +1] коэффициенты 
2 а b 
9| Ф-ла Гаусса 0,538469 0,906180 0,046910 0,2307655 
0,118463 0,239314 
10| Ф-ла типа | 97/180 = 163/180 = 17 183/360 
Гаусса = 0,53888 = 0,90555 1395 320805 
3799367 1345487 
11| To xe 6/11= 1011 = 1/22 5/22 
=0,545454 =—0,9090909 | 7381/64006= | 16577/69120— 
иу з ==0,1154- =0,24— 
+ 0,000328125 — 0,000170718 
12 › 0,55 0,91 0,0045 0,225 
0,11391936-- | 0,23910756 == 
је)“ =0,114— --0,239+ 


--0,00008064 + 0,00010756 


13| Ф-ла Шта- | 0,546377 0,9053237 0,04733815 0, 2268115 
2/17 4/17 


13" | Видоизме- | 0,5337205 0,9050822 0, 0474589 0,23313975 
нение ф-лы 3/25 6/25 
Штаермана 


15| Наивыгод- У 3/7= 
нейшая =, 
при 6=1 


*) См. сноску на стр. 40, 
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4* 


ординат Продолжениё 


Остаток 


[0; 1] и соответствующие 
при ординатах 


0,7692345 0,953090 
0,239314 0,118463 |0 


1/2 177/360 343/360 
71377321 | 0,238430397 | 0,118807949 | 0,00210 | 0,00490 | 0,0399 
249987721 _ 
=0,285523308 


1/2 17/22 21/22 

15643 4000-—= 0,239829282— 0,115328125=10,000355 |--0,00450| 0,0190 
—0'29— =0,24— | 0,1154 | 0,00271 |0,000480| 0,0133 
— 0,000314815 —0, 000170718 +0,000328126 


0,5 ‚115 0,955 | 
0,29394615— 10.23910756 = | 0,11391936— 0,00332 |-0.006285 0,0214 
294— „289+ [=0,114— |—0,0203 —0,00927] 0,0695 
0,00010756 — 0,00008064 


0,7731885 0,95266185 
0,235294 0,117647 0,0110 | 0,0214 | 0,0608 


0,9525411 | 
0,12 —0,00116 —0,00110] 0,0163 


| 
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Абсциссы на ООР 


20| Ф-ла типа 2/3 1 0 
Понселе **) 14/300 
(Ф-ла Харди) 
21 Кр. говая ИУ? 1 0 
а 5—1 / 0,5/15 


при них будут следующие: т 
у формулы 4“ при члене Са 


› › еле > съ 
3 

› > Це. И > съ 
> с" 

> > | A › ЕТ 
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*) У отмеченных этой сноской формул (см. прим. к табл. 2) в остат 


Пять 


Абциссы на отрезке 


коэффициенты 


Близкая 59/90 1 0 31/180 
к ней 1379/27714 43740000 
16078739 


0,17 
1/3,6877896 


0,18 
1/3,6274176 


0,175 
32/117,117 


1/6 
81 


козф. — 0,00159, а 
0,00312, 
0,00132, 

» — 0,100, 


> 


**) Получается из формулы Котеса на 7 ординат, если пропустить 


ординат Окончание 


Остаток 


[0;1] и соответствующие Ч 5 * 
Елон Код Алы ую козф аре коэф. при Уа 


ду | Ку | Мт 


12 149/180 1 
10443 — 0,27203626 | 0,04975824 | 0,00221 | —0,125 | —0,365 
=0,35641099 


0,5 0,83 1 
1178/3267— | 0,27116514 0,04854713 | 0,0131 | —0,1 —0,329 


==0,36057545 


0,5 0,82 1 
131/384= 0,27567821 0,05374887 | —0,0354 
=0,34114583 


0,5 0,8535534 1 
6/15=0,4 0,26666667 0,03333333 | 0,133 0,200 0,662 


ках появляются два дополнительных члена. Коэффициенты 


> 
при члене Ё ~ козф. — 0,00261 


О 

> ЗЕ а” 0,00440 
ГО 

> > ДА Др? 0,00209 
с 12 

> > Е > —0,06 


две ближайшие к середине ординаты. 
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бе“ 27: на отрезке 


1 | Ф-ла Чебышева | 0,266635 
0,422519 
1/6 


2 | Ф-ла Макло- 
рена . . . . 


16 36 56 


Що 


1280 


3| Круговая. . . арт о) ила 


=соз —— ш 


7 
0,62349 = 0,11346 
аз 
== 8 7 
4 | Ф-ла Гаусса . иа: Зар и 0,033765 


0,085662 


5 | Ф-ла типа Ta- 


усса....| 0,24 0,66 0,932 оди 


6| То же ....| 0,24 0,65 0,93 


0,24 0,66 0,94 ) 


Абсциссы на отрезке [0; 1] 


коэффициенты 


0,866247 |0,0668765 |0,2887405 |0,3666825 


1/6 1/6 


1/4 5/12 
247 254 
1280 1280 


0,188255 10,38874 


0,16340 |0,22314 


0,1693955 [0,3806905 
0,180381 [0,233957 


„17 0,38 
0,0864063 0, 1793362 02342575 


0,38 
77280 |0,2347141 


5 „38 
0,18558350/0, ‚23246287 
0,1855 0,2325 


0,033961 
1/12 


8 | Ф-ла Штаер- 


0,259565 0,932078 
0,681582 
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0,3702175 
3/12 


0,159209 
2/12 


Таблица 5 
ординат 


Остаток 


и соответствующие козф. при | козф. при 
при ординатах Аэ м ри 
Кэт пт! 


0,6333175 | 0,7112595 | 0,9331235 
1/6 1/6 1/6 


7/12 3/4 11/12 
0, 19843750] 0,19296875| 0, 10859375 | 0,143 0,369 0,757 


0,61126 |0,811745 | 0,950485 


0,22314 |0, 16340 0,11346 0,0156 0,0461 0,0884 


0,6193095 | 0,8306045 | 0,966235 
0,233957 |0,180381 | 0,085662 0 0 0 


0,62 0,83 0,966 
0,2342575 | 0,1793362 | 0,0864063 | —0,00191 | —0,00336 | —0,00376 


0,62 0,83 ga 
0,2347141 | 0,1777280 | 0,0875579 0,00125 0,00570 0,0128 


0,62 0,83 . [0,97 | 
0,23246287| 0,18553350] 0,08200363 | —0,0146 —0,0401 --0,0714 
0,2325 0,1855 0,0820 


0,6297825 |0,840791 | 0,966039 


3/12 2/12 1/12 0 0,00567 0,0133 
43 
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Шесть 


Е Абсциссы на отрез- Абсциссы на отрезке (0; 1) 
ју ке |— 1; Р коэффициенты 
2 


1 0 0,2 0,4 


9| Ф-ла Котеса .| 0,2 0,6 И 
19/288 75/288 50/288 


10| Наивыголней- 
шая при c=] | 0,285232 1 0 О, 1174724 |0,357384 


0,7650552 
0,0333335 |0,1892371 |0,2774294 
11 | Близкая к наи- 
выгоднейшей 
при е=И .|. 0,28 ‘0,276.44 0 0,12 0,36 
0,0341873 10,1922268 |0,2735859 


0 0,1175 0,3575 
„033306 |0,189392 |0,277302 


31 =» 277 16/21 1 0 5/42 5/14 
0,034623 |0,188536 10,276841 
| ла типа па 
| онселе . .| 1/4 3/4 1 0 3, 
| 103 22 ЕЧ 
3150 7105 225 
15 | Круговая при у 
с 1... .| 0,8090170:- 1 0 „0954965 |0,3454965 
==соз 72° 
"0,8090170 = 0,0200000 |0,1803715 |0,2996285 
== соз 36° 


Һір://гсіп.ого.р! 


ординат i Продолжение 


Остаток 
и соответствующие козф. при | козф. при | козф. при 
при ординатах а Г ће да 
т Ко МТ 
0,6 0,8 1 
50/288 75/288 19/288 —0,188 —0,520 --0,809 
0,642616 | 0,882527 1 
0,2774294| 0,1892371 0,0333335 0 0 0,0193 


0,64 0,88 1 
0,2735859] -0,1 0,0341873 0,00188 —0,0136 —0,0480 


0,6425 | 0,8825 1 
0,277302 | 0,189392 | 0,033306 0,000121 0,0003 0,0188 
9/14 37/42 1 

0,276841 | 0,188536 | 0,034623 | —0,00631 | —0,0176 | —0,0703 


0,2577778] 0, 0,0326984 0,0385 0,0576 0,0422 


0,6545035] 0,9045035 1 
0,2996285| 0,1803715] 0,0200000 0,00833 0,0438 0,0833 
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Семь 


k Абсциссм на отрезке Абсциссы на отрезке (0; 1] и 

Е 1; +1] коэффициенты 

2 а b с 

1 | Ф-ла Чебышева 0,323912 „0,883862 0,058069 | 0,2351715 | 0,338044 
сака и и и 


2 | Круговая.....| 0,3826834= 0,9238795= | 0,03906025 | 0,1464466 | 0,3086583 


= соз 450 0,0889822 | 0,1238095 | 0,1967322 


Ф-ла Гаусса .. 0,405845 Канео 2 0,025446 0,1292345 0,2970775 
0,064742 0,139853 0,190315 


4 |Ф-ла типа Гаусса 0 0,72 0,94 


0,03 0,14 0,30 
907457 236 | 0,13807608 | 0,18247433 
4° је) Тоже... . . 0,07: 0,13 0,1825 


КОЗУ, р а 0,0746 0,1381 ` | 0,1825 
5 ка 04 07 0, 0 3 
| “ ов ои Јава Ба | Мы 
в .....| од ол 05 | 0,026 0,13 од 
0,06421 0,14977 0,19083 
7 | Ф-ла Штаермана оома 0, 0,027777 0,1263015 0,297759 
115 215 3/15 
$ | Ф-ла Котеса... 13 23 1 0 1 1/8 
-41/840 9 9/280=» 
-216/840 | =27/840 
“| 1/20 1/4 
9 | Наивыгоднейшая 
при с=1..., 0,468850 1 0 0,0848885 | 0,265575 


10 | Близкая к наивы- 
годнейшей при . 
с. оо. во 0,47 0,83 0 
0,02397731 098803316 0,21547446 


*) См. сноску на стр. 48. 
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Таблица 6 
ординат 


соответствующие 
при ординатах 


0,5 0,661956 0,7648285 0,941931 

17 у 17 17 0,0260 0,0800 0,154 
0,5 0,6913417 0,8535534 0,96093975 
0,1809522 0,1967322 0,1238095 0,0889822 0, 00893 0,0277 0,0550 


0,7029225 
0,208280 0,190315 0,139853 


0,5 0,70 0,86 0,97 
020975446 | 018247433 | 0,13807608 | 0,07457236 | 0,000255 | 0,00441 0,0150 
0,2100 0,1825 0/1380 0,0745 0.00102 | 0,00312 0,00843 
0,2098 0,1825 0,1381 0,0746 |—0,000207 | 0;00186 0.00714 
0,5 0,7 0,875 0,975 
019717789 | 0:19946768 | 0/13971192 | 0,06223245 | 0,000548 | 0,01345 0,021 
0,5 07 0,87 0,975 
0,20438 0,19083 014277 0:06421 ,000288 |—0,0187 | —0,00349 
0,5 0,702241 
3/15 3/15 2/15 115 0,00412 | 0,0118 0,0216 
1 
мло ‚0593 „172 --0,321 
0: ,0765 |-0,196 —0,350 


0,9151115 1 
0,13841 


0,915 
0,13803316 


0,735 
0,21547446 


w vensen 


Работно тїї Natoma үшү ј а 
| 


ђје > 
WARSZA 


Bibli Мр: Мгсіп. -огд. pl- 


Маш» дан Абсциссм на отрезке (0; 1] 


b 


лизкая к наивы- 
годнейшей при 


*) При семи ординатах (см. прим. к табл. 2) в остатках формул, отме- 
енты при них будут соответственно следующие: 


hè № т 
ся ЕН “т 
У формулы 4’ 0,000600 0,00176 0,000832 
> > 4" —0,000246 —0,000286 —0,000315 
> > 8" —0,189 --0,0350 —0,0137 
> > . Ц"--0,000125 —0,000309 --0,00130 
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- 


ординат Продолжение 


1115 1112 1 
= | 0,21828531 = | 0,13893888= | 0,02264364 =| 0,00389 0,0116 0,0110 


168/705 | =0,21857183— = 0.18857143+ |--0,02285714- „000369 | 0,00109 | 0,00102 
+0,0002634= | --0,00029652 | +0,00036745 | —0,00021350 


0,14 0,92 1 
2-8 0,21918425 | 0,13383724 | 0,02173105 
-0,25049491 
56/900= | 0,21888889 | 0,13333333 | 0,02222222 |-0,00600 
—0,25111111 = 


1 13/14 1 
п 2206140 о 0,01794872 | 0,000583 0,0138 0,0333 
=0,26296293 


0,5 ү” 0,9 1 
дана 0,21849017 | 0,16189925 | 0,02809272 0,00567 |--0,00197 —0,0395 


‚183 0,218 0,1619 0,0281 0,00268 |—0,00115 —0,0200 
0,184 0,218 0,162 0,028 0,00421 0,000815 | --0,0175 


1 


0,2603175 0,2285714 0,1269841 0,0142857 0,00357 0,0179 0,0395 


ченных этой сноской, появляются три дополнительных члена; коэффици- 


h? hë h’ 
ст ЕБЕ бп 
У формулы 12 —0,000405 —0,00225 --0.00485 


14’ 0,000424 0,0000496 --0,0000766 
› > 14” 0,000640 0, 000640 —0,00184 
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LIAL 
ногАидоф 2 
онзиогинявйо 

еяоишо чизонго 


депнидеј 


0619969" 
1919690 


619679" 0 

1 169"0 
19999" 0 =2/% 
99619"0 
16519'0 
61189" 0 

0 

1 


гяонто 


06'0=49 66 0--2 


1202Ае.] епи чи-ф 


LILOJ ги-ф 
1=9 ер 
1=9 6/р=0 

98'0=9 %8'0=0 


10226 | сині м-ф 


еаәтмҹодә еи-ф 
22210} ви-ф 


6/1=2 
еззбе | пи ги-ф 
гнозпкио) ви-ф 


06/58 =v 
86'0=2 
9'0=2 


грозде | POHL чи-ф 


унпопеда ви-ф 


лениидо с 


чаениудо p 


маепиудо е 


чїрниїбо г 


вгениидо | 
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50000000 
2000000" 0 
110000000 "0 — 


9200000" 0 


0у60000“0 
84200000“0 
1100000“0 
8610000“0-- 
#%00000 "0— 
9602000 ‘0— 


6970000“0 

е600000"0-- 
е00000“0 
810000“0 


22210 
ноиќидоф 5 
оняиздинявдо 
ечоито *LHIOHLO 


еядишо 


с619969“0 
87619969“0 
19:619969"0 


сссд9во“0 


8219929" 0) 
сессд9е9“ 0) 
Р1с9929*0 
66689290 
6199290 
1999690 


99999290 
с0199с9*0 
1759989*0 
99899890 


[=> 
[=> 
уб'0=> 


[> 

1-2 
$6‘0=2 
та 


9/6=9 
'8'0=4 


она 
:229Ае.1 епн чи-ф 


"01/49 ) 


“02 
КА 
зени\ о L 


222104 тиАидоф 


ује=4 


'91'0=4 85'0=2 
'99'0=4 %с'0=р 


NÁR епиі чи-ф 


евопичооћ Ри-ф 
29910} ви-ф 


14 
1=9 
16:04 
6:04 


1 
| 
УП => 
99:04 фа“? | тениидо 9 


осо 
090 
бо"? 


1ениидо $ 


әинәжготой || 
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Таблица Та 


Расположение формул в порядке уменьшения 
даваемой ими абсолютной ошибки для у == ѕіп х 


Ф-ла трапеций (2 ординаты) 
> одной ординаты 
› типа Гаусса с абсциссами 
› Симпсона (3 ординаты) 
› типа Гаусса с абсциссами 
> > > > > 
» Котеса (4 ординаты) 
» типа Гаусса с абсциссами 
> > > > > 
› Котеса (5 ординат) 
> Чебышева (4 ординаты) 
» Котеса (6 ординат) 
» типа Гаусса с абсциссами 


> > > > > 
> > > > > 
> > » > > 
> > › › › 
› › › › › 
› › › › > 


> Котеса (7 ординат) 

› Типа Гаусса с абсциссами 
> > > > > 

» Чебышева (6 ординат) 

> типа Гаусса с абсциссами 


а=0,0 
а=0,6 (2 ординаты) 


а=0,58 (2 ординаты) 
а=52/90 (2 ординаты) 


а=0,5; b=1 (4 ординаты) 
а=7/9 (3 ординаты) 


а=4/9; b=1 (4 ординаты) 
а=0,34; b=0,86 (4 ординаты) 
а=1/4; 6=3/4; с=1 (6 ординат) 
а=2/3; b=1 (5 ординат) 
а=0,24: ђ=0,66; с=0,94 

(6 ординат) 
а=0,5; b=0,9 (5 ординат) 
а=0,65; 6=1 (5 ординат) 


а=0,55; 5=0,91 (5 ординат) 
а=0,52; #=0,9 (5 ординат) 


а=0,28; 6=0,76; с=1 (6 opan- 
нат) 


а=0,24; 6-0,66; с=0,93 
(6 ординат) 


а=0,4; $=0,8; с=1 (7 орди- 
нат) 

а=7/15; ђ=5/6; с=1 (7 opan- 
нат) 


а=0,4; 6 0,72; с=0,94 (7 op- 
динат) 
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взоншо изоно 


зиножиогод | | 


ї=2 19/9=49 ср г? 
„Гас? 8:04 +'0=2 
Р6'0==2 :51'0==9 усо 


тезоће | епил чи-ф 
22210} еиќидоф 


1=9 „Зад ФР 
1=2 1910 =9 '85'0=2 
$6'0=2 :99'0=4 #20 => 
760--2 '99'0==@ с 0-- г 


182242 | ециј чи-ф 


газат ги-ф 


54 


ст.огд.р 


И 


http 


Таблица Ва 


Расположение формул в порядке нъшения даваемой имя 
не форма ошибки Ала у == 1--соз 2х 


Ф-ла трапеций (2 ординаты) 
› одной ординаты а=0,0 
› Симпсона (3 ординаты) 
> типа Гаусса с абсциссами 
> > > > > 
> › > >+ > 
И Котеса (4 ординаты) 
› типа Гаусса с абсциссами 
> › › > > 
> > > > > 
› Котеса (5 ординат) 
› Чебышева (4 ординаты) 
> Котеса (6 ординат) 
› типа Гаусса с абсциссами 


> > > > > 
> > > › › 
› › › › › 
> > > > › 
› › › › › 


> > » > >» 

» Котеса (7 ординат) 

› Чебышева (6 ординат) 

» типа Гаусса с абсциссами 
> 

> 

> 


а= 0,6 (2 ординаты) 
4--0,58 (2 ординаты) 
a= 52/90 (2 ординаты) 


а= 0,5; 6 =1 (4 ординаты) 
а==7/9 (3 ординаты) 
а == 4/9, 6 == 1 (4 ординаты) 


4--0,28; 6--0,76; с=1 
(6 ординат) 

а--2/3; 6-- 1 (5 ординат) 

а= 0,24; 2 — 0,66; с = 0,94 
(6 ординат) 

а = 1/4; b=3/4; с=1 (6 op- 
динат) 

а= 0,24; 5—=0,66; c = 0,93 
(6 ординат) 


а--0,34; 5=0,86 (4 орди- 
наты) 


а=0,5; b= 0,9 (5 ординат) 


http: И 


а=0,65; 6=1 (5 ординат) 
з 2" "а 4--0,55; 6 =0,91 (5 ординат) 
2 АТ. 4 а= 0,52; b= 0,9 (5 ординат) 
з ди та > а=0,4 20,8; с=1 (7 op- 
динат) 
z + га а а=7/15; b= 5/6; с =1 (7 op- 
динат) 
» , ад э а= 0,4; 6 = 0,72; САХ 
(7 ординат) 
| 
5* 55 
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Таблица 9 


5 Ошибка 
Черт. ih ~ , 
1 орди- а=0,0 0,576 
ната 
Ф-ла трапецић —0,424 
ото | Ф-ла типа Гаусса —0,191 
а == 52/90 
Ф-ла Симпсона 0,242 
3 Span | Ф-ла типа Гаусса 0,00579 
а =7/9 
Ф-ла Котеса 0,0627 
Ф-ла типа Чебышева 
а = 0,16; 6 = 0,8 0,0134 
4 Ф-лы типа Гаусса: 
opan- |} 4--31/90; b= 78/90 —0,00535 
и а == 61880; b = 155/180 0,00422 
а = 1/3; b 53790 0,00270 
а= 45: b= „0569 
а = 4/9; уз: --0,0510 
Ф-ла Котеса —0,0397 
5 орди- Ф-лы типа Гаусса: 
нат а = 47/90; b = 81/90 0,000319 
а = 97/180; b = 163/180 0,000612 
Ф-ла Котеса —0,0199 
Ф-лы типа Гаусса: 
в} 0,00164 
6 е b= 0, , 
но ои с= 0.94 
a= 0.24 
b = 0,66 0,000019 
с=0,93 
Ф-ла Котеса 0,00544 
Ф-лы типа Гаусса: 
60:72 | 0,000025 
- , , 
ра“ c= 0,94 : 
а= 7/15 
b= 59 0,000158 
гуаг 


Таблица да 


Расположение формул в порядке уменьшения даваемой ими 
абсолютной ошибки для у == sin? х 


Ф-ла одной ординаты 
трапеций (2 ординаты) 
Симпсона (3 ординаты) 
типа Гаусса с абсциссами а==52/90 (2 ординаты} 
Котеса (4 ординаты) 
Котеса (5 ординат) 


Котеса (6 ординат) 

типа Чебышева с абсциссами а=0,16; b= 0,8 (4 ординаты) 
› Гаусса с абсциссами а=0,45; 6—1 (4 ординаты) 
› > > > : a= 4/9; b= 1 (4 ординаты) 
> > тъ › a= 1/9; (3 ординаты) 

Котеса (7 ординат) 

типа Гаусса с абсциссами а = 31/90; а” (4 орди- 

наты 


> > 6 > а = 61/180; b = 155/180 (4 op- 
динаты) 


a= 1/3; b = 78/90 (4 ординаты) 


а--0,24; 6--0,66; с=0,94 
(6 ординат) 


а ==97/180; b = 163/180 (5 op- 
динат) 


а = 47/90; 6 = 81/90 (5 орди- 
нат) 

a= 1/15; 6=5/6; с =1 (7 op- 
динат) 


a= 0,4, b= 0,72; с = 0,94 
(7 ординат) 

a= 0,24; b= 0,66; с = 0,93 
(6 ординат) 


а Три, уу 
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Таблица 10а 


Расположение формул в порядке уменьшения даваемой ими 


абсолютной ошибки для у = 


ла трапеций (2 ординаты) 
одной ординаты 
Симпсона (3 орлинаты) 
Котеса (4 ординаты) 
типа Гаусса с абсциссами 
Котеса (5 ординат) 
типа Гаусса с абсциссами 
Котеса (6 ординат) 
типа Гаусса с абсциссами 
Котеса (7 ординат) 
типа Гаусса с абсциссами 
> > › › 
› › › › 
Чебышева (5 орлинат) 
типа Гаусса с абсциссами 


> > > > 


Чебышева (4 орлинаты) 
типа Гаусса с абсциссами 


Чебышева (6 ординат) 
типа Гаусса с абсциссами 
> > › 
› › › 
› › › 


› 


1—х? 


а=0,0 


а ==4/9; В==1 (4 ординаты) 
a= 52/90 (2 ординаты) 
а=0,65; # = 1 (5 ординат) 


а=2/3; b= 1 (5 ординат) 
а=0,5 6—1 (4 ординаты) 
а=0,6 (2 ординаты) 


a= 0,28: 0,16; е==1 
(6 ординат) 
а =17/9 (3 ординаты) 


а = 1/4; 6В= 8/4; с=1 
(6 ординат) 
а=0,4; 5=0,8; е =! 
(7 ординат) 
а =7/15; b= 5/6; с=1 
(7 орлинат) 
а=0,34 b= 0,86 (4 opan- 
наты) 


а= 0,52; b= 0,9 (5 ординат) 
а= 0,5; 6 =0,9 (5 орлинат) 
а == 0,58; b= 0,91 (5 ординат) 
а = 0,24; b= 0,66; с = 0,93 
(6 ординат) 
а=0,4 b= 0.72; с = 0,94 
(7 ординат) 


a= 0,24; b= 0,66; е = 0,94 
(6 ординат) 
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Таблица 11 


Черт, 13. 


1 ордината а=0,0 


Ф-ла трапецић 

Ф-ла Гаусса, Чебышева 

Ф-ла та Гаусса 
4--0, 


2 ординаты | 
| 


Ф-ла Симпсона 
Ф-ла Чебышева 
Ф-ла Гаусса 
Ф-ла типа Гаусса 


Ф-ла Котеса 
Ф-ла Чебышева 
Ф-ла типа Гаусса 


3 ординатъг 


4 ординаты |< Ф-ла Гаусса 
а=0,34; 0--0,86 


Ф-ла Котеса 

Ф-ла Чебышева 

Ф-ла Гаусса 

Ф-ла типа Гаусса 
а=0,52; b= 0,9 


5 ординат 


Ф-ла Котеса 

Ф-ла Чебышева 

Ф-ла Гаусса 

Ф-ла типа Гаусса 
2) 

Ф-ла Котеса 

Ф-ла Гаусса 

Ф-ла типа Гаусса 
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aA Е Таблица Па 


Расположение формул в порядке уменьшения даваемой ими 
абсолютной ошибки для 
а Утв 


Ф-ла трапеций 
_ Симпсона 
Котеса 
одной ординаты а =0,0 
Гаусса (2 ординаты) 
типа Гаусса с абсциссами а=0,6 (2 ординаты) 
Чебышева (3 ординаты) 
Гаусса (3 ординаты) 
типа Гаусса с абсциссами а 7/9 (3 ординаты) 
Чебышева (4 ординаты) 
Чебышева (5 ординат) 
типа Гаусса с абсциссами a = 0,34; b=0,86 (4 ординаты) 
Гаусса (4 ординаты) 
Чебышева (6 ординат) 
типа Гаусса с абсциссами а= 0,52, b= 0,9 (5 ординат) 
Гаусса (5 ординат) 
Гаусса (6 ординат) 
типа Гаусса с абсциссами г. а=0,24 b = 0,74; е =0,95 
(7 ординат) 
Гаусса (7 ординат) 


типа Гаусса с абсциссами а=0,24; 0--0,66; с = 0,94 
(6 ординат) *) 


*) Этот результат, будучи лучшим, чем даёт точная формула 
Гаусса на 7 ординат, является случайным для данного частного 
примера. 
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Таблица 12 


Черт. 14, 


Ошибка 


1 орди- ГО : 
ната а==0,0 0,015625 0,001819 
Ф-ла трапецић 0 —0,138 
2 ор- Ф-лы типа "у аусса: 
динаты а=0, 0,016640 0,00283 
а=0,58 0,016971 0,00317 
3 ор- Ф-ла Симпсона 0,0104167 |—0,00339 
динаты Ф-ла ен Г Гаусса : 0,013033 |—0,000773 
Ф-ла ана 0,013095  —0,000771 
4 ор- > лы типа 05 
динаты = 0,34; 0 0,86 0,013485 | —0,000321 
20:45; та _-- | 0,01495 0,00115 
Ф-ла Котеса 0,014713 0,000907 
5 ор- Ф-лы типа Гаусса: : і А 
динат а= 0,52 2=0,9 0,0137796 |—0,0000262 
а=0,5 6=0,9 0,01381698 | 0.0000111 
а = 59/90; b=1 0,013949 0,000144 
Ф-ла Котеса у 0,0143815 | 0,000576 
Ф-лы типа Гаусса: 
69 а=0,24 
> b= 0,66 0,013790697; —0,0000151: 
| е=0,93 
а =0,24 
b= 0.66 0,0138369 0,0000311 
c= 0,94 | 
Ф-ла Котеса 0,013946 0,000140 
| Ф-лы типа Гаусса: 
7 ор- а=0,4 0,01379745 
динат 5=0,726 - : РА 
= 0,94 
а =7/15 
= е 0,0138177 0,0000119 
је: 


"һїр;//гсїп.огд.р! 


Таблица 12а 
Расположение формул в порядке уменьшения даваемой ими 
абсолютной ошибки для y= |x Их (1 — х)]° 


Ф-ла трапеций (2 ординаты) 
> (Симпсона (3 ординаты) 
» типа Гаусса с абсциссами 
> > пр 45 
> одной ординаты 
» типа Гаусса с абсциссами 
> Котеса (5 ординат) 
» типа Гаусса с абсциссами 
» Котеса (4 ординаты) 
> Котеса (6 ординат) 


» типа Гаусса с абсциссами 


› › › > › 
> Котеса (7 ордннат) 


› типа Гаусса с абсциссами 


> > › > > 
› › › › › 
› › › › › 
› > > > > 
› . " » ? 


а=0,58 (2 ординатм) 

a= 0,6 (2 ординаты) 
а=0,0 

a= 0,45; 6 = 1 (4 ординаты) 


а =7/9 (3 ординаты) 


а=0,34; 8--0,86 (4 орди- 
наты) 


а = 59/90; b= 1 (5 ординат) 


a= 0,24; b= 0,66; с = 0,94 
(6 ординат) 


a= 0,52; b= 0,9 (5 ординат) 


a = 0,24; 6 = 0.65; с = 0,93 
(6 ординат) 


ае; 026/6 с==1 
(7 ординат) 


а—0,5; b= 0,9 (5 ординат) 


а==0,4; b= 0,72; с=0,94 
(7 ординат) 


һїр;//їсїп.огд.р! 


Таблица 13 
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8 Ошибка 
Черт, 15- 
| 
1 орди- а=0,0 1,407195 —0,946 
ната 
Ф-ла трапеций 5,5 3,147 
2 орди- Ф-лы типа Гаусса: 
наты а== 0,6 3,0929 0,740 
в =0,58 2,0099 —0,343 
Ф-ла Симпсона 2,6048 0,252 
Ф-ла пеций, повторённая 
$ орли- два Жы а 3,4536 1,101 
57" Ф-ла типа Гаусса 
а= 7/9 2,2912 —0,0615 
Ф-ла Котеса 2,5407 0,188 
Ф-ла трапеций, повторённая“ 
три раза 2,86948 0,517 
У Ф-ла Чебышева 2,31425 —0,0385 
В ја Ф-лъ типа Гаусса: 
a= 61/180; b= 155/180 2,33797 —0,0147 
а=1/3; 678/90 13/15 | 2,39232 0,0396 
а = 4/9; в-- 1 2,34963 —0,00308 
5 Ф-ла Котеса 2,33086 —0,0219 
SpA Е: Ф-ла трапеций, повторённая 
четыре раза 2,6322 0,280 
65 


1 


Ф-ла Симпсона, повторени 
два раза \| 

Ф-лы типа Гаусса: 

a= 47/90; 6-- 81/90 

a= 97/180; b = 163/180 


Ф-ла Котеса 


Ф-ла трапеций, повтореннал 
пять раз 


Ф-ла Гаусса 
Ф-лы типа Гаусса | 


Ф-ла Котеса 


Ф-ла трапеций, повторённая 
шесть раз 


Ф-ла Симпсона, повторённая 
три раза 
Ф-лы типа Гаусса: 


„Продолжение 


2,358404 


2,38412 
_2,36970 


2,311443 


2,492406 
2,3679295 


2,41284 


` 2,276406 


2.289401 
2,44506 


2,30804 


http://rcin.org.pl 


Ошибка 


. 0,00569 


0,0314 
0,0170 


—0,0413 


0,1397 
0,0152 


0,0601 
--0,0763 
—0,0633 


—0,0924 


—0,0447 


Таблица 13а 


Расположенне формул в порядке уменьшения даваемой ими 
абсолютной ошибки для ~ 


Ф-ла трапеций (2 ординаты) 


> 
> 
> 
> 
> 
> 


трапеций, повторённая два раза 

типа Гаусса с абсциссами 

4 > > > 

трапеций, повторённая три раза 

типа Гаусса с абсциссами 

трапеций, повторённая четыре 
раза 

Симпсона (3 ординаты) 

Котеса (4 ординаты) 

трапеций, повторённая пять раз 

трапеций, повторённая шесть раз 

типа Гаусса с абсциссами 


Котеса (7 ординат) 
типа Гаусса с абсциссами 
> > > > 


Симпсона, повторенная три раза 
Котеса (6 ординат) 
типа Гаусса с абсциссами 


Чебышева (4 ординаты) 
типа Гаусса с абсциссами 


> > > > 
Котеса (5 ординат) 
типа Гаусса с абсциссами 


> > > > 


Гаусса (6 ординат) 
типа Гаусса с абсциссами 


> > > > 


Симпсона, 
раза #) 


повторвнная два 


(8 ординаты) 
а--0,0 (1 ордината) 
а= 0,6 (2 ординаты) 

(4 ординаты) 

а = 0,58 (2 ординаты) 


(5 ординат) 


(6 ординат) 
(Т ординат) 
a= 0,24; b= 0,66; с = 0,93 
(6: ординат) 


а =7/9 (3 ординаты) 
а= 0,24; b= 0,66; с = 0,94 
(6 ординат) 
(7 ординат) 


а = 1/3; аза (4 орли- 
наты 


а = 47/90, 6-- 81/90 (5 op- 
динат) 
а ==4/9; 6—1 (4 ординаты) 


а = 97/180; b= 163/180 
(5 ординат) 

а =7/15; В--5/6; с=1 
(7 ординат) 


а — 61/180; b= 155/180 
(4 ординаты) 

а==7/15; 0 5/6; с=1 
(7 ординат) 


(6 ординат) 


*) Этот результат является шуй для данного частного 
_ случая. 


http:// 


rcin.org.pl 


учар“ Ун Р 


Таблица 14 


И 
49 
as зедбобвоща 
47 
4 
44 8 
4 
42 
y 
Черт, 16 
В Руде а--0,0 0,710634 
Ф-ла трапецић 0,55 
раче Ф-лы типа Гаусса: 
жи а=0,6 0,63744 
а = 0,58 0,64204 
| Ф-ла Симпсона 0,657089 
Ф-ла трапеций, повторённая 
3 орди- 2 раза 0,630317 
наты “а, типа Гаусса 
а << 1/9 0,645424 
Ф-ла Коте 
че Дре че повторвнная 0546186 
? раза у 0,64231 
Ф-ла Чебышева 
4 орав Ф-лы типа Гаусса: 0.9991 
а = 61/180; b= 155/180 0,646330 
а= 1/3; Б = 78/90 == 13/15 0,647342 
а= [9761 0,649492 
$ Ф-ла Котеса 0,649053 
д“ Ф-ла трапеций, повторённая 
четыре раза 0,64474 


Ошибка 


Продолжение 


8 Ошибка 


а=0,4 
а 
с=0,94 


р = 5/6 


: 0,6467906 | —0,0000109 
ф== 


Ф-ла Симпсона, `повторённая 
два раза 0,649555 0,00275 
5 орди- Ф-лы типа Гаусса: 
да a= 47/90; -- 81/90 0,646588 |--0,000213 
| ta= 97/180; b = 163/180 0,646640 |—0,000161 
Ф-ла Котеса 0,6484018 | 0,00160 
Ф-ла трапеций, повторённая 
пять раз 0,64587 |—0,00093 
Ф-лы типа Гаусса: 
а=0,24 
повин Ер 0,646899 | 0,000098 
с=0,94 
а=0,24 
b= 0,66 0,646757 |—0,000044 
с=0,93 
Ф-ла Котеса ` 0,647566 0,000755 
Ф-ла трапеций, повторённая 
рт раз 0,646681 0,000120 
Ф-ла Симпсона, повторённая 
три раза 0,648139 0,00133 
Ф-лы типа Гаусса: 
"7 opan- 
нат 0,6467849 | —0,0000166 


Ба 
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Таблица 14а 


Расположение формул в порядке уменьшения даваемой ими 
абсолютной ошибки для Y= Âp 


Ф-ла трапеций (2 ординаты) 


одной ординаты 
трапеций, повторённая 
раза 

Симпсона (3 ординаты) 
типа Гаусса с абсциссами 
Котеса (4 ординаты) 
типа Гаусса с абсциссами 
трапеций, повторённая три 
раза 

Симпсона, повторённая два раза 
типа Гаусса с абсциссами 
Котеса (5 ординат) 

трапеций, повторённая четыре 
раза 

Котеса (6 ординат) 

типа Гаусса с абсциссами 
Симпсона, повторённая 3 раза 
Чебышева (4 ординаты) 
трапеций, повторённая пять раз 
Котеса (7 ординат) 

типа Гаусса с абсциссами 


два 


трапеций, повторённая шесть раз 
типа Гаусса с абсциссами 


а=0,0 (1 ордината) 


(3 ординаты) 
а==0,6 (2 ординаты) 


a= 0,58 (2 ординаты) 


(4 ординаты) 
(5 ординат) 
а = 4/9; b= 1 (4 ординаты) 


(5 ординат) 


а= 7/9 (3 ординаты) 
(7 ординат) 


(6 ординат) 


а = 1/3; 6 = 78/90 
(4 ординаты) 

а = 61/180; 5 = е 
(4 ординаты 


а = 97/180; 6 = 163/180 
(5 ординат) 
(7 ординат) 
а=0,24; b= 0,66; е =0,94 
(6 ординат) 
а= 0,24; b= 0,66; с = 0,93 
(6 ординат) 
а= 0,4; 6 = 0,72; с = 0,94 
(7 ординат) 
а= 1/15; b = 5/6; с = 1 
(7 ординат) . 
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egom 
-чдәһ 


ст впимоет 


160 =2 
соста 
КАЛЕ 


-дов ту јижеззитодез 


(езодв jeun) 
гнияниътод 


"же Кидоф оп кеннәгоныче “ëg 
тенипдо у-ви 


+4 
с 


езо/е | 


ЧУЕ Фор Се оре 
сес90*0 | 86790'0 | 21190:0 
ере | 90:94 680 
"46100 | 9ев/0“0 | 195/0'0 
Фе | Живи 909102 
181600 | бвебо“0 | $7860'0 
Фо ОТ чај 261 
96110 |сипо| 116110 
Форс 408:9 ЕИ) 
оссго свето 91.0 
%9'1 %® %8'© 
988/1'0 | 829/1°0 | 8990/10 
Фото | %80'1 НАЧИ 
100000 | 28612:0 | 61660 
%71 ‘0 %ус0 | 61:0 


же 06| = гядишо 
рес60“0 


9000,0 
08“ еб| == вядишо 


9080,0 |. зебого 
ос '01| = езбито 


10860 "0 с0211'0 
906 01| = гядишо 
шо севето 
Фос с | = гядито 
сар] 0 свет" 0 
9091 | = ехдито 
698/1'0 58181 ‘0 
%0 == РЯ0ИИЮ 
cazzo | 55555'0 
%9 == Рядишо 


11282" 0 


евош Да ще енең | ©0=гр 

Е е Жым :ийводиподео| vIe dig 

ао рол ока "ЧЕНИ (взе решиш) | зиноьене 
найее вто сенки | оонвој 


мегкидоф оп веннэионыча “Ро 
лениџдо с=и 


== 6 енвидэе: п] Фэки || 


d7, 
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